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Abstract 

The interference of polariton guiding modes propagating through "traffic circle" waveguides 
composed of dielectric spherical particles is investigated. The dependence of intensity of the wave 
on the position of the particle was studied using the multisphere the Mie scattering formalism. 
We show that if the frequency of light belongs to the passband of the waveguide, electromagnetic 
waves may be considered as two optical beams running along a circle in opposite directions and 
interfering with each other. Indeed, the obtained intensity behavior can be represented as a simple 
superposition of two waves propagating around a circle in opposite directions. The applications of 
this interference are discussed. 

PACS numbers: 
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I. INTRODUCTION 



Low-dimensional aggregates of nanoparticles absorbing, scattering and guiding electro- 
magnetic waves are attracting growing interest because they can be used in various micro- 
and nano-systems where optical ei ierg y is received, transferred and converted on a subwave- 



length scale 
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12| . If the material's refractive index is sufficiently 



large, then these aggregates can possess long-lived quasistates, which are polariton modes 
representing the superposition of light and polarization oscillation of the material 
These modes propagate within the aggregates, while their losses due to light emission are 
very small. Consequently, low- dimensional aggregates of particles can be used as waveguides 
for optical energy transport. This transferred energy can stimulate a number of processes 
such as the photoexcitation of a molecule or a chemical reaction 3, Q] • 

For example, if we excite the polariton mode within a linear chain of particles using 
a point harmonic source located near the end of the chain, then there exists a frequency 
domain, i. e. the passband, where almost all of the energy from the source will be guided 
by the chain [l^ . Other modes in the related frequency domain can be treated as bound or 
guiding modes. To our knowledge, the first nano-waveguide capable of transferring optical 
energy a distance of 100 nm was constructed by Atwater and coworkers linear 
chain of spherical silver or gold particles. The use of metallic particles has an advantage in 
that the losses due to photon emission can be almost completely suppressed |l2|, yet the 
efficiency of energy transfer is suppressed due to optical energy absorption by free electrons. 

Guiding systems of dielectric particles have negligible absorption losses, while radiative 
losses can also be avoided as dielectric materials have characteristically larger refractive 
indices with ra^ > 2 j^, [l^. The finite size of the guiding energy band in periodic arrays of 
particles results in the the formation of slow light modes at the band edge 12(]. Therefore, it 
ma y b e more convenient to use dielectric materials having very weak light absorption [o], Q, 
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12| . Frequency pass bands of linear chains were investigated H 
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l^and it was 



found that the top of the frequency pass band corresponds to the top of the Brillouin band 
edge in quasi-momentum space, while the bottom of the frequency pass band is determined 
by the guiding mode criterion (see Eq. (fT^ below) j^, 1^. It was discovered that the closer 
the frequency of electromagnetic waves emitted by the source to the upper boundary of the 
pass band and the longer the chain, the more effective the waveguide and the lower the 
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amount energy is dissipated into the surroundings. It was shown that for +00, the 

rate of decay of electromagnetic waves decreases as 7(A^) ~ A^~^ 

Another example of a low- dimensional aggregate of particles possessing wave-guiding 
properties is the circular array of spherical particles. Because of its symmetry, one can 
expect that the highest quality factor can be attained for particles arranged in a circle. A 
circle has no sharp ends in contrast with a particle chain where the lifespan of polariton is 
limited to the time it takes to travel from one end of the chain to the other [lO^. It has been 
demonstrated that the bound modes in a circular array of spherical particles possess a 



quality factor that grows exponential 
arrays of cylinder- shaped antennas ( 



ly wit 
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1 the number of particles, just as in the circular 

,0,Q). 



As previously mentioned, waveguides comprised of spherical particles can be used to 
manipulate optical energy on a subwavelength scale. In addition, modes at the edge of the 
passband possess a vanishing group velocity which creates an interesting opportunity to 
realize the phenomenon of slow light [12]. The properties of polariton modes within pass 
bands are quite similar to those of electrons within a wire. The question of interest is whether 
we can exploit several properties known about electron behavior in a wire for wave-guiding 
systems. In this manuscript, we investigate the effect of interference on the transport of 
optical energy. The geometry of the waveguide shown in Figs. [T] [2], [3] is similar to the one 
illustrating the Aharonov-Bohm effect in electronic systems. The sample geometry can be 
described as a "traffic circle" without an entrance or exit (Fig. [1]) or with them (Figs. EllHl). 

The polariton is excited by a point source located near one of the aggregate particles as 
shown in Figs. [H El [3l Practically, this source can be created by placing a dye molecule 
within an end particle (particle A) py]. When the frequency of the source belongs to the 
passband, almost all of the optical energy generated by the source is transferred through the 
first chain |l^. After traveling through the first chain, an electromagnetic wave enters the 
circle at the junction between the circle and the first chain (particle B). It then splits into two 
beams, one of them propagating around the circle clockwise and the other counterclockwise. 
These two beams rejoin at the junction between the circle and the second chain (particle C) 
and travel through the second chain to particle D. 

Actually, the beam propagating along the circle clockwise from particle B to particle C 
may pass through particle C from one segment of the circle to the other segment and continue 
traveling from particle C back to particle B. Analogously, the other beam propagating along 
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the circle counterclockwise from particle B to particle C also may pass through particle C 
and continue traveling back to particle B in the other segment of the circle. This recircula- 
tion in both the clockwise and counterclockwise directions remarkably enhance interference 
effects. Therefore, we investigate the interference of waves having frequencies corresponding 
to resonances of circular arrays where this enhancement is maximized. 

To describe our expectations about the interference effect, we can suggest the following 
simple scenario perfectly that is also applicable to an electronic system. Each of two beams of 
light may be expressed as a wave exp{ikx} if the beam travels clockwise and ex]i{ik{L — x)} 
if the beam travels counterclockwise. Here, k is the wave vector of the beam and it is defined 
by the frequency of the source, x is the coordinate of a point on the circle (we suppose that 
for point B x = 0) and L is the circumference of the circle. So, if these two beams interfere, 
the intensity of the field is given by the superposition of these two waves: 



exp{ikx} + exp{ik{L — x)}. 



(1) 
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FIG. 1: Traffic circle without chains composed of spherical particles. 
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This interference may have a number of interesting apphcations including, for instance, fil- 
tering and selective enhancement of a guiding signal. It is not clear though whether this 
simple interference model Eq. ([T]) is applicable to optical modes because of the strong spatial 
dispersion and essential non-local coupling of spherical particles composing the waveguide. 
In this manuscript, we show that despite the above-mentioned problems, our simple inter- 
ference model Eq. ([1]) works quite well. 

To ascertain that interference of polaritons in a circular array takes place, we have pro- 
vided numerical investigations of a circle with and without attached chains. We have studied 
three particular circle without chains, a circle with chains situated diametrically 

opposed to each other and a circle with chains situated orthogonally to each other. For each 
of these three cases, we have investigated two situations: the presence and the absence of 
absorption. 

The system under consideration is a traffic circle composed of Nc = 80 identical dielectric 
spherical particles (see Figs. [H El [3l). The radii of particles were chosen to be unity. The 
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FIG. 2: Traffic circle with collinear chains composed of spherical particles 
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chains (linear arrays externally joined to the circle) are composed of same particles. The 
number of particles in these chains is Ni = N2 = 20. The angle between the chains is 
denoted a. If the chains are collinear, then one has a = Fig. [21 and if the chains are 
orthogonal to each other then one has a = 7r/2 Fig. [31 

We introduce an oscillating magnetic dipole oriented orthogonally to the system plane as 
the source of electromagnetic waves. A magnetic dipole is chosen because it corresponds to 
;he guiding band with the lowest frequency and, consequently, the best guiding properties 
9j. Indeed, as shown in the mentioned work, the lowest frequency guiding mode originates 
from magnetic-like Mie resonance. If the system under consideration is a circle with chains, 
the source is situated near the end of one of the chains, as shown in Figs. [21 [3 If the circle 
without attached chains is considered Fig. [H the source is situated outside of the circle 
near the first particle. In this manuscript, we consider GaAs particles having a refractive 
index = 3.5. GaAs is a convenient to consider since it possesses one of the highest 
refractive indices amongst conventional optical materials and, consequently, it has better 




FIG. 3: Traffic circle with perpendicular chains composed of spherical particles 
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guiding properties (cf. [9!]). The absorption effect was also probed and we discuss it when 
reporting the results of our calculations. 

All calculations were performed using the multisphere Mie scattering formalism ( [lol. I20I ) . 
This formalism is introduced in Sec. [ITl where we also discuss the calculations of the polariton 
dispersion needed to determine the resonant modes of the traffic circle. This formalism is 
used in Sec. IIIII to investigate the interference of polaritons in different geometries. We 
investigate the intensity of the polariton around the circle in arrays shown in Figs. [H El [3] 
as well as the intensity transferred through the waveguide with chains at arbitrary angles. 
The obtained dependence was compared with Eq. ([T]) representing the sum of two optical 
beams to verify if the dependence may be approximated with this function. 

The dependence of the intensity of an electromagnetic wave scattered by the a spherical 
particle in the circle on an arbitrary point on the surface of this particle is demonstrated 
in Figs, m O [6], [TJ [HI It was found that in all investigated cases the dependence may be 
approximated by Eq. 

In Sec. llVl the results are summarized and their applications to various systems of interest 
are discussed. 

II. MULTISPHERE MIE SCATTERING FORMALISM. 
A. General formalism 
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The multisphere Mie scattering formalism [19|, |20|] (see also 
been developed to study the scattering of electromagnetic waves in aggregates of spherical 
dielectric particles. The electromagnetic field outside each particle may be represented as 
the sum of all incident waves falling on this particle (say, the j-th particle) and of diverging 
waves scattered by this particle: E{r) = El^^{r) + Ei^^{r), H{r) = Hl^^{r) + Hla{r). 
The incident wave Ej^^{r), Hl^Jj") may be expressed as a linear combination of spherical 
vector harmonics Nml{r), Mmi{r) taken with respect to the center of j^^ particle 

+00 n 

Eincir) = Y.Y. ' (PlnN^lir - r,) + qLM^l{r - r,)) , (2) 

n=l m=—n 
+00 n 

Hincir) = E E {I'mnN^lir - r,) + p^^^M^lir - r,)) . (3) 

n=l m=—n 
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Analogously, the scattered wave Ei^^{r), Hi^^{r) is a linear combination of spherical vector 
harmonics Nrnn{r), Mmn{r) corresponding to outgoing waves 



hoo n 



EUr) = E E ^ KnNl^nir - r,) + V^M^l{r - r,)) , (4) 

71=1 m=—n 
+00 71 

HUr) = E E KnN^^nir - r,) + a^^Mlir - r,)) . (5) 

n=l ?Ti=— 71 

Here is the coordinate of center of the j^^ particle; indices j = 1,...N enumerates all 
particles, N is the total number of particles in the array; index n = l,...,+oo represents 
angular momenta of spherical vector harmonics, and index m = —n,...,n stands for the 
projection of angular momentum on the z-axis. The coefficients a^„, p^Ti^ Qmn 
partial amplitudes of scattered and incoming spherical waves at the j*'* sphere. 

The incident wave falling on the j-th particle is the sum of the external waves emitted by 
the source and the waves scattered by all other particles. Because of the boundary conditions 
at each sphere, amplitudes of these waves obey the relationship 

N +00 1/ 
N +00 u 

where a„ and bn are Mie scattering coefficients depending on the size and the refraction index 
of particles. They originate from boundary conditions for an electric field between particles 
and air. Parameters A'-^,^^^ and -B^n/xi^ vector translation coefficients that depend on the 
relative position of spheres j and /. These coefficients express delayed multipole interactions 



of different spheres in the frequency domain 



9|. These translation coefficients can be used 



to represent the outgoing wave scattered by the l^^ particle as a sum of spherical vector 
harmonics Nml{r), Mml{r) taken with respect to the j*^ particle as 

+ 00 71 

Ml'Jir - n) = E E {A'in,uM^l{r - r,) + Bli^^,N^l{r - r,)) , (7) 

71=1 m=—n 
+00 n 

N^ir - n) = E E {Blin,uMil{r - r,) + A'i^^,N^l{r - r,)) . (8) 



71=1 m=—n 



To complete the setup of the problem Eq. ([6]), one has to determine the coefficients p^Ti 
and pi^j^ for the incoming wave of the point magnetic dipole source having the coordinate Vq 



8 



(see Figs. [H El [3]). This external wave can be expressed as a linear combination of spherical 
vector harmonics Nmn{r), Mmnir) 

+00 n 

EeMr) = E E * {Pn.nN^l{r - ro) + Q^nM^f) (r - ro)) , (9) 

n=l m=—n 
+00 n 

H,.Ur) = E E {Q-^nN^li^ - ro) + P^„M(fl(r - ro)) . (10) 

n=l m=—n 

Here Pmn and Qmn are partial amplitudes of an emitted wave characterizing the point 
source. As mentioned above, we chose the source to be the magnetic dipole parallel to the 
z direction and perpendicular to the x — y sample plane. Then the wave emitted by the 
source can be described by the parameters Pmn and Qmn defined as Pmn = Qmn = for 
n> 1, alsoP_ii = Poi = -Pii = 0, Q-ii = Qii = 0, and Qoi = 1. 

The coefficients p4,n and are associated with coefficients Pmn and Qmn by means of 
the relations 

+00 u 

Pmn ~ ^ ] ^ ^ {A'mnfj.uP^i.i' + ^mufiuQ fJ.)^) ' (H) 
u=l ^=—v 
+00 V 

Qmn ~ ^ ] ^ ^ {■^mnfiu^fJ-i' ~^ ^mn^iuQ f^v) 5 (12) 

U=l fl= — l' 

where vector translation coefficients are taken between the source and the center of the 
corresponding sphere. Thus, using Eqs. (fTT!) and (fT2|) one can calculate the right hand side 
of Eq.®. 

Eq.(l6]) has the form of a system of linear equations, with partial amplitudes and 
tPmn yet to be determined. Our main task is to calculate these partial amplitudes and to 
investigate their dependence on the number of spheres. 

In this work, the dipolar approximation is used which is sufficiently accurate to charac- 
terize the system 



ihe dip 



12| . Using this approximation, Eq.(l6]) can be expressed as 

<ii _L V A^^ -\- R^^ U \ - R^^ 

y '=T^'^T ^ ^ ^. (13) 

The numerical solution of these equations is used later in Sec. Illll to investigate interference 
of waves in the circular traffic waveguide Figs. [H [2], [31 
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B. Determination of the Dispersion Law 



Partial amplitudes and characterize interference of two polaritons running clock- 
wise and counterclockwise. The parameter k entering into Eq.([T]) is the quasi- wave vector 
of a polariton and it is determined by the frequency u of the wave emitted by the source. 
The relationship between u and k is the wave dispersion law which has to be determined 
for the analysis of our simulation data. Below we briefly described how we did that. 

If the radius of the circle is large enough as in our case (a circle containing = 80 
particles), one can use the dispersion law for an infinite linear chain. Any frequency lj 
lying in the conducting band of an infinite linear chain corresponds to a certain value of the 
quasi- wave vector k. 

The frequency domain of bound modes, or the pass band, is the domain of frequencies 
where losses due to photon emission are completely suppressed in the limit N ^ oo. This 
frequency domain exists where photon emission is forbidden by the momentum conservation 
law. To illustrate, consider a polariton within an infinite chain of particles parallel to the z- 
axis. According to the light dispersion law, we can define the wave vector of light in vacuum 
as g = ujc, where c is the speed of light. If an optical mode can exist outside of the chain as 
a free photon it still has the fixed projection of the wave vector to the z-axis given by k so 
we should have Qz = k due to translational symmetry. Then the perpendicular component 
of the electric field wave vector that describes the photon emission outside of the system 
can be expressed as q± = \/q'^ — k"^. For unbound modes, this projection is real, meaning 
that the optical energy can escape from the chain of particles, while for bound modes it 
is imaginary so that the mode intensity decays at distance p from the chain as e~^v^^^~?. 
Thus modes with uj/c <k are bound, while the modes with uj/c> k are unbound. Besides 
that, according to the Bloch theorem the absolute value of k cannot exceed the value vr/a, 
where a is the period of the chain. Therefore, bound modes can possess quasi-wave vectors 
k from the interval 

^<k<^-. (14) 

c a 

Certainly, this condition may be satisfied only if uo/c < ix/a. This guiding mode criterion 
is equivalent to the condition that the interparticle distance a should be less than half the 
resonant wavelength: a < X/2 [2^. 

To find the dispersion law for a polariton within an infinite chain , l5| one has 
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to calculate the eigenvalues of the frequency uo for the given quasi- wave vector k. The 
eigenfrequencies can be computed using the condition of the existence of nontrivial solutions 
of a homogeneous system of equations 

N +00 u 

N +00 u \^'^) 
binn _|_ ( Alj U _|_ Ulj \ = {) 

According to the Bloch's theorem, the partial amplitude in an infinite periodic chain can be 
represented by 

= a p*''^ h> = h e'^"" (16) 

mn ^mn'^ ; "mn "mn'^ ; \ 

where k is the quasi- wave vector and x = aj is the coordinate of the center of j*'* particle 
in the chain (remember, we set a = 1). 

Substituting Eq. ([16]) into Eq. ([15]), we get (cf. Refs. Q,^) 

+ 00 

~^ 5^ (^?Tmmi/(^; ^) ^^mu ~l~ -B^nmui^^ j ^) ^mu) 0) 

'Too (17) 

^^^T^ ~l~ i^^mnmui^^ y ^) ^mi/ ~l~ BYnnmui}^ i ^) ^mu) 0- 
i/=l 

We used the momentum projection m conservation law due to the rotational invariance of 
the system. 

In the dipolar approximation, Eq. fll7l) takes the form 

+ 00 

^TJT + XI (^mlml('^; ^) O-ml + B„ilml{^ i k) h„il) = 0, 



^^=1 (18) 



If m = ±1, the modes are called t-modes (meaning transverse modes with polarization 
perpendicular to the chain), and if m = the modes are called /-modes (meaning longitudinal 
modes with polarization parallel to the chain). In this manuscript we are interested in t- 
modes. 

As was discussed above, the source of electromagnetic waves in our problem of interest is 
a magnetic dipole directed perpendicular to the plane of the circle. This source excites only 

transverse electromagnetic waves. Therefore we are interested in equations for m = 1 

+00 

^+Y: (Anii(^, k) an + 5nii(^, k) 6n) = 0, 



+00 

61 M 



(19) 



^ + E (^1111 (c^, k) fell + Bnnioo, k) an) = 0. 

^ ' u=l 
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This system contains only two variables an and bn and the dispersion law is defined by the 
following condition 

k)] (-^ + Auuiuj, k)] - {Buuiu;, k)f = 0. (20) 



Eq. (l20l) can be solved only numerically. We used the Newton- Raphson method to find 
quasiwavevector for certain frequencies of interest following Ref. Qj. 

III. RESULTS 

In this work three arrays of particles were considered, including the circle of particles 
without attached chains Fig. [1] and two structures with chains Figs. [2] and [3l Each of these 
structures is investigated for particles of ideal non-absorbing material. Our preliminary study 



of wave transmission through a traffic circle at various frequencies has demonstratec 



that it is 
These 



most efficient for modes corresponding to resonances of a circular array (see Ref. {9]). 
resonances takes place at frequencies corresponding to integer projections q = —N/2, N/2 
of a polariton's angular momentum to the z-axis. These modes can be characterized by 
quantized quasi-wavevectors k = 2'Kq/ L. The upper boundary of the pass band corresponds 
to g = N/2, i.e. k = it. For the circular array under consideration this upper boundary is 
given hj uj = 0.8434. 

In this work we investigate a circular array composed by = 80 particles. To characterize 
the interference, we chose two modes having angular momenta qi = 38 and q2 = 37 and 
frequencies Ui = 0.8424 and 0^2 = 0.8413, respectively, located at the top of the passband 
{q = 40). Interference behavior is quite similar for other modes with q > 30. The mode at 
the top of the band characterized by g = 40 does not show any interference because the wave 
intensity does not change from sphere to sphere; though, the field changes its sign between 
neighboring particles. 

The intensity of a polariton mode can be characterized by squared partial amplitudes 
a{^^ and bl^^. Since we are using the dipolar approximation, only coefficients a{^^ and fo^n 
for n = 1 should be considered. The partial amplitudes a^^, b''^^ and are all equal to 
zero because of the planar structure of our system. Therefore, the intensity of a polariton 
within the j-th particle is characterized by three coefficients a^n, an and b^^. For the sake 
of simplicity, we define it as Ij =\ p. Using the multisphere Mie scattering formalism 
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|l9l |. we investigated the dependence of intensity of the wave on the j-th particle on the 
position of the particle. The position of the particle is defined by its number j and may be 
described as the length of the segment of the circle between the first particle and the j-th 
particle counted counterclockwise: Xj = {j — 1). 

The intensity of the wave should be compared with the one following from Eq. ([T]) 

Ij oc| ex.p{ikx} + exp{iA;(L — x)} |^ /4 = cos^(7rg — kx) = cos^ {27i qx/N). (21) 

Since for guiding modes of interest one has q ~ N/2, this function strongly oscillates for 
non-integer x and is difficult to show in our graphs. Therefore, we compare it directly to 
our simulations in Fig. H] only, while in all other graphs it is replaced with the function 

Ij oc cos2(27r(iV/2 - q)x/N), (22) 

which coincides with Eq. (pT]) for all integer points x of interest. For chosen modes with 
g = 38 and q = 37, this function can be expressed as cos^(2 ■ 27rx/40) and cos^(3 ■ 27ra;/40), 
respectively. 

Below we report our calculations for the electromagnetic field within dielectric spheres for 
two selected modes for all three configurations shown in Figs. [H [21 [3] and then investigate 
the effect of interference on the transmitted intensity. 

A. Circle without attached chains 

Consider the intensity behavior in a circle without attached chains (Fig. [1]). In Fig. H] 
and Fig. O we show the results of numerical calculations of the wave intensity Ij depending 
on the number of particles for two different resonant source frequencies Ui and uj2 defined 
previously. Solid curves on Figs. HJ [5] represent the data fit by functions Eqs. ( 12T|) and ( l22l) . 
respectively, with a numerical prefactor chosen as a scaling factor to make the fit intensity 
identical to the calculated one for j = I. According to Figs. [Hand [5] at both frequencies ui 
and UJ2, there is a perfect agreement between our calculations and our simple interference 
model Eqs. ([21]), Thus guiding modes in a circular array excited by a point source can 

be treated using the interference of two optical beams propagating along a circle in opposite 
directions. 

We verified that the weak absorption Im Ur < 0.01 does not qualitatively affect the in- 
terference behavior Eqs. ([22l) . but remarkably reduces the wave intensity. This is because 
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guiding modes of circular arrays possess extremely high quality factors [95 growing expo- 
nentially with the number of particles. Even absorption as weak as Im rir ~ 10~^ results in 
remarkable broadening of all resonances. 



B. Traffic circle with attached chains 



Consider the traffic circle with attached chains Figs. [21 [3l The circle is composed of 
Nc = 80 particles, and the chains are composed of Ni = N2 = 20 particles. The source 
of electromagnetic waves is again a magnetic dipole oscillating normally to the plane of the 
circle. It is located close to the end of one of the chains (see Fig. [2], [3]). 

As in the previous case, we investigate the dependence of intensity of a polariton at 
particle j on the number j. The influence of the chains on this dependence is examined. We 
begin our consideration with the system with the chains situated coUinearly as in Fig. [3 

The results of calculation are presented in Fig. [6l Symbols denote the intensity Ij 
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■ Fit of intensity by IG)=5.69-10'cos'(2k(41-j)), k=38-27i/80 
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FIG. 4: High resolution representation of intensity of the wave depending on particle number j 
within a circle without attached chains 
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calculated numerically, and solid lines show fit using the analytical function Eq. (122|) with 
arguments Xj = 2(j — 1). The squares and the upper solid line correspond to the frequency 
ui = 0.8424 and the triangles and the lower solid line correspond to the frequency lj2 = 
0.8413. 

It is clear that the symbols are located near the solid lines. Therefore, the dependence 
of intensity Ij j can be approximately expressed by the function Eq. (!22l) . similarly to the 
traffic circle without attached chains. However, the agreement of the numerically calculated 
intensity with the analytical expression Eq. (122|) is less accurate than in the case of the 
circle without attached chains. This is clear because guiding mode propagation in a circle 
with chains is much more complicated than in the circle without chains due to scattering at 
its junctions. 

Consider a circle with chains forming right angles with each other. In Fig. [TJ the results of 
numerical calculations of the absolute value of the intensity Ij are shown for the two modes 
under consideration. In this case, similarly to the case of collinear chains, the dependence of 
intensity on the position of a particle in the circle may be well approximated by the function 




FIG. 5: Intensity of a wave depending on particle number j within a circle without attached chains 
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Eq. (122|) . Therefore, the approximation of intensity by the superposition of two waves is 
acceptable for any position of the second chain. 

The intensities for a circle without chains and for a circle with chains in the absence of 
absorption differ by six orders of magnitude. This effect is due to the reduction of the mode 
quality factor due to polariton scattering by the chains. Indeed, quality factors of resonant 
modes of circular arrays increase exponentially with the number of particles in the circle. 
However, quality factors of circles with chains behave similarly to those of finite linear chains 
where they increase with the number of particles according to the power law Six orders 
of magnitude difference in intensity is also due to a difference in quality factors of resonant 
modes under consideration. 

Another interesting phenomenon is that absorption in a circle with chains has a weaker 
influence on intensity. This is another consequence of the difference in quality factors. 




FIG. 6: Intensity Ij versus particle number j within a circle with collinear chains 
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C. Transmitted intensity affected by interference 



Traffic circle interference is characterized by its transmitted intensity. This is the intensity 
at the end of the second chain (particle D - see Figs. [2], [3]). We investigate its dependence 
on junction position, which can be any particle in the circular array. Similarly to Figs. 
[21 [3] the second chain is placed in a radial direction perpendicular to the circle, while its 
starting point can any arbitrary particle within the circle. This dependence characterizes 
the interference effect on the output signal. 

We expect that the dependence can still be approximated (at least roughly) by the func- 
tion Eq. 0221) assuming that the intensity at the end of the second chain should be propor- 
tional to the intensity of field at its origin. Since the intensity in the j-th particle is defined 
by the interference of polaritons and may be approximated by Eq. fl22l) . one can assume 
that the intensity at the end of the chain behaves similarly. 

In Fig. [8], the calculated dependence of the transmitted intensity Ij on the junction 
position j is shown. As shown previously, the intensities calculated numerically are marked 
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by symbols and analytical fits by Eq. (122|) are shown using solid lines. One can see that 
the dependence of partial amplitude 601 at the end of the second chain on the junction 
position can be approximated by the function Eq. (1221) . Yet this approximation is rough. 
For several positions of the junction the error can be as high as 30%. Still this approximation 
is convenient for a qualitative description of behavior of an electromagnetic field at the end 
of the second chain. The infiuence of absorption on the transmitted intensity of field at the 
end of the second chain is weak similarly to a circle without chains as considered previously. 

IV. CONCLUSION 

In this work we have investigated the interference of guiding modes in a traffic circle 
waveguide of dielectric spherical particles using the multisphere Mie scattering formalism 
at wave frequencies corresponding to two resonant guiding modes. For different sample 
geometries we demonstrated that the field distribution within the circle as well as the trans- 
mitted intensity may be well approximated by interference of two optical beams propagating 
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clockwise and counterclockwise around the circle. 

The main conclusion of this work is that strong interference exists in the waveguides 
composed of particles and is quite similar to interfenrence as seen in electronic systems. 
This interference can be used to control the output (transmitted) intensity through varying 
the system geometry and/or the mode frequency. These results must also be applicable 
to the subwavelength waveguides consisting of metal particles 3, S] despite absorption by 
conducting electrons and can be used to narrow the transferring signal in the frequency 
domain. 

It is difficult to apply our results directly to experiments with whispering gallery modes 
in arrays of polystyrene particles {l| because the dipolar approach is irrelevant for these 
modes having large angular momenta. However, interference effects might take place there 
because coupling of sphe.e= i= e^entially local Another interesting application can be 
developed for the propagation of microwaves in arrays of coupled antennas [3|, [ISj which 
are very similar to our particle aggregates. Interference can be used in this application for 
frequency selective emission, absorption or transport of the microwave signal. 

We do not discuss in detail frequency dependence of the transmitted intensity. Our 
preliminary study shows that this dependence is very complicated because the interference 
of clockwise and counter clockwise beams is superimposed with circular array resonances. 
This problem is still under investigation and we do not have a reasonable interpretation for 
it yet. 

This work is supported by the Air Force Office of Scientific Research (Award no. FA 9550- 
06-1-0110). Authors acknowledge fruitful discussions with Vasilii Asratov, Olga Samoylova, 
Mark Sulkes and Arthur Yahjian. 



[1] V. N. Astratov, J.P. Franchak, S.P. Ashili, Appl. Opt. 85 (23), 5508 (2004) 
[2] L.I. Deych and O. Roslyak, Phys. Rev. E 73, 036606 (2006) 

[3] R. A. Shore and A. D. Yaghjian, "Traveling electromagnetic waves on linear periodic arrays 

oflossless spheres". Electron. Lett. 41, 578-580 (2005). 
[4] S. Fan, J. N. Winn, A. Devenyi, J. C. Chen, R. D. Meade and J. D. Joannopoulos, "Guided 

and defect modes in periodic dielectric waveguides," J. Opt. Soc. Am. B 12, pp. 1267-72 



19 



(1995). 

[5] R. D. Meade, A. M. Rappe, K. D. Brommer, J. D. Joannopoulos, O. L. Alherhand, "Accurate 
theoretical analysis of photonic band-gap materials", Phys. Rev. B 48, 8434 (1993). 

[6] Z. Y. Tang, N. A. Kotov, "Onc-dimcnsional assemblies of nanoparticles: Preparation, prop- 
erties, and promise," Adv. Mater. 17, 951-962 (2005). 

[7] S. A. Mayer, P. G. Kik, H. A. Atwater, S. Meltzer, E. Harel, B. E. Koel, A. A. G. Requicha, 
"Local detection of electromagnetic energy transport below the diffraction limit in metal 
nanoparticle plasmon waveguides," Nature Mater. 2, 229-232 (2003). 

[8] S. A.Mayer,M. L. Brongersma, P. G. Kik, S.Meltzer, A. A. G. Requicha, B. E. Kocl, H. A. 
Atwater, "Plasmonics A route to nanoscale optical devices," Adv. Mater. 15, 562-562 (2003). 

[9] A. L. Burin, "Bound whispering gallery modes in circular arrays of dielectric spherical parti- 
cles," Phys. Rev. E 73, 066614 (2006). 
[10] A. L. Burin, G. C. Schatz, H. Cao, M. A. Ratner, "High quality optical modes in low- 
dimensional arrays of nanoparticles. Application to random lasers," J. Opt. Soc. Am. B 21, 
121-131 (2004). 

[11] M.I. Gozman, I.Ya. Polishchuk, A.L. Burin, Optical Modes in Linear Arrays of Dielectric 
Spherical Particles: A Numerical Investigation. Proceedings of the 9th International Confer- 
ence on Transparent Optical Networks, v. 4 136 - 139 (2007) 

[12] Gail S. Blaustein, Michael I. Gozman, Olga Samoylova, I. Ya. Polishchuk, and Alexander L. 
Burin, Guiding optical modes in chains of dielectric particles. Optics Express, Vol. 15, Issue 
25, pp. 17380-17391 (2007) 

[13] Michael I. Gozman, I. Ya. Polishchuk, and Alexander L. Burin. Features of Propagation of 
Light in the Linear Array of Dielectric Spheres Proceedings of the 10th International Confer- 
ence on Transparent Optical Networks, (accepted for publication) (2008) 

[14] Michael I. Gozman, I. Ya. Polishchuk, and Alexander L. Burin Light Propagation in Linear 
Arrays of Spherical Particles. Phys. Let. A (in press) (2008) 

[15] R. W. P. King, G. J. Fikioris, R. B. Mack, Cylindrical Antennas and Arrays, Cambridge 
University Press, 2005. 

[16] D. K. Freeman, T. T. Wu, IEEE Trans. Ant. Prop. 43, 340 (1995). 

[17] V. M. Apalkov, M. E. Raikh, B. Shapiro, Phys. Rev. Lett. 92, 253902 (2004). 

[18] E. I. Smotrova, A. I. Nosich, Opt. Quant. Electronics 36, 213 (2004); S. V. Boriskina, 



physics/0509 203, submitted to IEEE J. QUant. El. (2005). 
[19] Y. L. Xu, "Electromagnetic scattering by an aggregate of spheres: far field," Appl. Opt. 36, 
9496-9508 (1997). 

[20] Y. L. Xu, "Scattering Mueller matrix of an ensemble of variously shaped small particles," J. 
Opt. Soc. Am. A 20, 2093-2105 (2003). 

[21] H. W. Ehrespeck, H. Poehler, "A new method for obtaining maximum gain from Yagi anten- 
nas," IEEE Trans. Antennas Propag. AP-7, 379-386 (1959). 



21 



